Computational Investigation of Experimental Interaction Impedartce 
Obtained by Perturbation for Helical Traveling-Wave Tube Structures 

Carol L. Kory /// - 

ANALEX Corporation/NASA Lewis Research Center ^ 

21000 Brookpark Road, MS 54-5 J? *7 V / 

Cleveland, Ohio 44135 

James A. Dayton, Jr. 

NASA Lewis Research Center 
Cleveland, Ohio 44135 

ABSTRACT- Conventional methods used to measure the cold-test interaction impedance of 
helical slow- wave structures involve perturbing a helical circuit with a cylindrical dielectric rod 
placed on the central axis of the circuit. It has been shown that the difference in resonant 
frequency or axial phase shift between the perturbed and unperturbed circuits can be related to 
the interaction impedance. However, because of the complex configuration of the helical circuit, 
deriving this relationship involves -several approximations. With the advent of accurate three- 
dimensional helical circuit models [1, 2, and 3], these standard approximations can be fully 
investigated. This paper addresses the most prominent approximations made in the analysis for 
measured interaction impedance by Lagerstrom [4] and investigates their accuracy using the 
three-dimensional simulation code MAFIA. It is shown that a more accurate value of interaction 
impedance can be obtained by using three-dimensional computational methods rather than 
performing costly and time consuming experimental cold-test measurements. 

I. INTRODUCTION 

Obtaining slow-wave circuit interaction impedance is an important part of the design 
process for a traveling-wave tube (TWT) as this parameter is related to the gain and efficiency of 
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the tube or the minimum length of a backward-wave oscillator. On-axis interaction impedance is 
proportional to the strength of coupling between the RF wave and the electron beam, defined for 
the n' 1 ' RF space harmonic [5] as 
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where E^O) is the on-axis longitudinal electric field magnitude of the n th axial space harmonic, P 
is the time averaged RF power flow and p„ is the axial phase constant of the n* space harmonic 
defined by 
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where P is the fundamental axial phase constant and L is the pitch length. 

Synchronism with the beam for the electromagnetic wave propagating on the helix is 
desirable for the fundamental RF space harmonic; thus, we need to calculate the fundamental (n 
= 0) harmonic term in equation 1 to determine the on-axis interaction impedance which is related 
to the gain of the circuit and is used as an input in various TWT modeling codes such as [6]. To 
determine the interaction impedance of the backward wave mode which is related to backward 
gain (backward wave oscillations) we need to calculate the impedance associated with the n = -1 
space harmonic is required. 

Typically, interaction impedance is obtained by performing cold-tests on the helical 


circuit or a scale model of the circuit. This experimental procedure is time-consuming and 
expensive, but until recently it was thought to be the most accurate means of determining 
interaction impedance. Most often, measurements on helical circuits are performed by using 
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nonresonant perturbation techniques where the axial phase shift is measured at a certain 
excitation frequency for a perturbed and an unperturbed circuit. After some lengthy algebraic 
manipulation, Lagerstrom shows that this difference in phase shift, A|3, is proportional to the 
perturbed and unperturbed electric fields integrated over the perturbing rod volume divided by 
the power within the structure (equation 3. 1 8 of [4]) 

£ 0 £ o (e'- 1) j*E .E* dv 

A/3 = ^ . 3 
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E* denotes the complex conjugate of the original electric field and E' denotes the perturbed field. 
If representations of the unperturbed and perturbed fields are known, equation 3 can be solved 
providing an expression relating this measurable quantity to the interaction impedance. 

However, because of the complexity of the helical structure, several approximations become 
necessary in representing E* and E\ It is well known that these assumptions will cause some 
error in the measurement, and several researchers have addressed this issue in an attempt to 
improve upon Lagerstrom’s work [7], [8]. Nevertheless, several approximations are still 
involved making direct calculation of equation 1 with a three-dimensional computer code [1,2, 
and 3] the most accurate method. 

An accurate, three-dimensional, computational helical cold-test model has been 
developed and reported on in [1]. This model allows, for the first time, a complete investigation 
of the conventional approximations which go into helical circuit analysis. A summary will be 
given of Lagerstrom’s analysis emphasizing the major assumptions made, followed by the results 
from an investigation of each of these assumptions using MAFIA. MAFIA (Solution of 
MAxwell’s equations by the Finite-Integration-Algorithm) is a powerful, electrodynamic code 
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that is used for computer-aided design of two-dimensional and fully three-dimensional 
electromagnetic devices [9, 10]. 

In addition to the approximations that go into the experimental impedance calculations, 
there are several human factors which can cause inaccurate measurements. For example, the 
perturbing rod may be displaced from the central axis of the circuit slightly, or the dielectric 
constant of the perturbing rod may be uncertain. Also, reflections may be present because of 
mismatches at the input/output couplers, nonuniformity in the helical pitch, or because of other 
manufacturing errors. 

The circuit used as a model for this study is the helix slow- wave circuit from a 40 Watt, 
18-40 GHz TWT for the millimeter-wave power module (MMPM) [11]. The experimental 
circuit includes a rectangular, tungsten, helical tape supported by T-shaped BeO rods inside a 
conducting barrel as shown for the end view in Figure 1. The specific results presented here 
apply only to the mentioned circuit; however, these results can serve as a general guide for 
similar devices, and the computational techniques are readily applicable to other TWT’s. 
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Figure 1 Hughes MMPM helical slow-wave circuit 

II. THEORY OF INTERACTION IMPEDANCE MEASUREMENTS [4J 

Lagerstrom uses the perturbation method to distort original fields in the circuit in order 
to gain knowledge about the strength of the fields in the region of interaction. Reviewing the 
perturbation formula of equation 3, we recognize that it is necessary to derive expressions for 
unperturbed and perturbed electric fields within the perturbation rod volume. If these 
expressions can be derived in terms of the magnitude of the n lh axial space harmonic of the on- 
axis longitudinal electric field, E^O), the interaction impedance can be expressed in terms of the 
measurable quantity Ap. Thus, as a basis for the helical field configuration we start by 
considering a circular aperture with radius a which identically repeats itself every distance L in 
the axial direction. The space harmonic nature of the field solutions is such that any field 
quantity can be written as 
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where the time variation e i “‘ is understood and \ is the n* Fourier amplitude coefficient which is 
a function of r and 0 only. Inside a homogeneous region of a circular aperture in free space, the 
space harmonic components are solutions of the homogeneous wave equation in cylindrical 
coordinates. So, the longitudinal electric field can be expressed as 
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E z (r,0,z) = e-^E z (r,0,z) = e~^ + 5 
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where I v (y n r) and K^,(y n r) are the modified Bessel functions of the first and second kind and y n is 
the radial propagation constant of the n ,h space harmonic defined by 

Vn = \P\ • 1(2 > k2=0)2 ^o 6 

where co is the radian frequency 2nf, and p 0 and s 0 are the free space permeability and 
permittivity, respectively. When the axis is included in the circular aperture, K^r) cannot 
represent a physical field, forcing b TO to equal zero. Thus, in the region 0 < r < a 
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The phase velocity v pM of the vn ,h space harmonic is given by v pW =— . The interaction 

Pn 


impedance for the v« ,h space harmonic at a non-zero value o: ‘ radius can be defined as 
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where P denotes the time-aveiaged real power flowing in the mode 
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of which this vn ,h space harmonic is a part. Since there can be no azimuthal variation on the axis 
of a circular aperture, the interaction impedance at r = 0 reduces to equation 1 . 

Although at any position all of the space harmonics are present and contributing to the 
total field, the perturbing rod can be strategically placed in a region where the space harmonic of 
interest is dominant. It is obvious that on the axis of a circular aperture, the 0-varying space 
harmonics vanish (v = 0) and the n = 0 space harmonic will dominate. If a dielectric rod of 
infinitesimal diameter were placed on the axis, only this space harmonic would be present in the 
region of perturbation. Accordingly, only this space harmonic would be involved in the 
perturbation integral of equation 3 making the problem straightforward. An infinitesimally thick 
rod is, of course, not physically possible, and in practice the rod must be thick enough to produce 
a measurable perturbation. Therefore, other space harmonics will be present over the rod volume 
making the problem far more complicated. In [4], the perturbing rod is assumed to be small 
enough in diameter that the space harmonics other than v=n= 0 can be neglected over the entire 
perturbation volume. Another assumption made is that the TM portion of the fields dominates 
the TE portion. Therefore, electric fields can be expressed in terms of the axial and radial 
components only, neglecting the azimuthal component. 

The approach used in [4] to calculate the local fields for the perturbation integral 
considers that the perturbation rod is placed in a region containing a “snap-shot” of the original 
field configuration. Lagerstrom refers to this as the “incident” field. The “induced” or 
“reflected” field, AE, is calculated as the additional field needed to satisfy the perturbing rod 
boundary conditions. With these calculated, E' can be calculated for use in equation 3. 
Lagerstrom justifies this procedure as follows: the propagation constant is changed only slightly 
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by the addition of the dielectric rod, and since the media of the slow-wave circuit determines the 
unperturbed propagation constant, the fields at the perturbat on would be expected to change 
only slightly as well. 

Applying Lagerstrom’s approximations that only the v=«=0 space harmonic will be 
present and the TE portion of fields can be neglected within the perturbation volume, and the 
perturbed and unperturbed fields are equal outside the perturbation rod, we proceed in solving 
equation 3. To express the fields in the actual structure, [4] begins as if to solve the boundary 
value problem of a perturbing rod of radius b in a sheath helix of average radius a and axial turn 
length L. Since the particular configuration of the circuit is not involved, the local field solutions 
at the rod and the evaluated integral should apply for an arbitrary structure with a dominant non- 
azimuthally varying harmonic. 

The TE and TM fields can be expressed separately in terms of the longitudinal fields, Ez 

A 

and Hz. Using the notation E = e~^ E where E is now a function of r only. Maxwell’s 
equations [12] are used to give 
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From equation 7, the unperturbed longitudinal electric field in the region inside the helix, 0 < r < 
a, can be expressed as 

E z =AI Q { Y r) 9 
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Since / 0 (y r ) = 1 on the axis, A is the magnitude of the axial electric field on the 
axis,£ z (0) =A. From equation 8, the transverse fields needed in the perturbation integral of 
equation 3 associated with the TM portion of the fields can be derived: 

P 
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Next, a cylindrical dielectric rod of radius b and dielectric constant e'e 0 is placed on the 
helical axis, consequently creating a new boundary condition. Two regions now exist, I, the 
region between the rod and the helical conductor and II, the region within the dielectric rod (See 
Figure 2). 



Figure 2 Helical perturbation rod model 

Using equations 5 and 8, the fields in the two regions can be expressed as 
Region I (b < r < a): 
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Region II (0 < r < b): 
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For the perturbed case, the magnitude of the longitudinal field on the axis is 

Erf/iO) = A n = 2T' Z (0) . At the rod surface, the normal D field and tangential E field must be 

continuous, so 

E a (b) =E ai (b ) 

E r l(b) = S E r /j (b) * 13 

As mentioned previously, the rod is assumed to be placed ir a snap-shot of the unperturbed 

t i 

fields, so the incident fields are equal to the original fields, y * y ,and P~P at the perturbation 

rod surface. With these assumptions, E zI can be extended to the axis and evaluated resulting in A 
= A, = E z (0). The perturbed on-axis longitudinal field magnitude, E' z (0) = Ezii(O) - A,, 
Substituting equations 1 1 and 12 into 13 and eliminating B I; Lagerstrom derives the correction 
factor p, relating the perturbed and unperturbed on-axis longitudinal electric field magnitudes as 
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Substituting the electric field components and taking advantage of the assumption that 
there is no azimuthal or axial variation, equation 3 can be expressed as 


u 

2kcos 0 {£'-\) z o(r)E z o(r) + E r o(r)E* T o(r)^dr 


A/J = 


4 P 


18 


The fields within the integral can be expressed in terms of the on-axis fields using equations 9, 
10 and 12: 
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Lagerstrom defines two additional correction factors: 
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To simplify the notation p e is defined as 


, p 2 . i Vo*) 

Pe - P\iP2+ ' Pv - ' ,1 r 
y y s - 1 _ , _ 


r ,+ • > 
u eV -V 


L-fi) 

eF -F 


22 


and the integral of equation 1 8 becomes 
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Substituting equation 23 into equation 3, Lagerstrom arrives at 
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Consequently, the on-axis interaction impedance can be expressed in terms of the change in 
propagation constant as 
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It should be noted, as Wang and Carter [7] point out, that Lagerstrom omitted a term 
when evaluating p e (equation 3.27 of [4]). They report in [7] that this omission results in an error 


12 



which is larger than the measurement error at higher frequencies. Equation 22 gives the correct 
formula. 


The above analysis which applies to the nonresonant perturbation method can also be 
applied to the resonant perturbation method where the resonant frequencies are measured in a 
truncated section of circuit before and after perturbation. The change in frequency can be related 
to the change in axial phase shift by equation 3.54 of [4] 

26 
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where p, v p and f are the unperturbed axial phase shift, fundamental phase velocity and 
frequency, respectively, and v g is the group velocity of the perturbed structure. 

Another common approximation made when the dielectric rod is small, is that the axial 
electric field is uniform in the transverse direction and the transverse electric fields can be 
neglected. Making these assumptions simplifies the problem to a great extent since the field 
tangential to the perturbing rod must be continuous through the boundary; thus, the magnitude of 
the perturbed field is equal to that of the unperturbed field everywhere [4]. In this case, equation 
3 is easily solved and the interaction impedance is expressed as 

= £ z (0)£;(0) _ 2 A/7 

0 2p 2 P coP(e'-X)e 0 nb 2 P 

Lagerstrom refers to equation 26 as the first order perturbation expression. He has improved 
upon this theory by deriving the correction factor p e to account for the transverse fields, the 
variation of the axial fields in the transverse direction and the difference between the perturbed 


13 



and unperturbed fields. This first order perturbation expression is still used, however, by various 
researchers so it is worthwhile comparing these results with the corrected theory by Lagerstrom. 


III. SIMULATED INTERACTION IMPEDANCE 

The on-axis interaction impedance for the n th RF space harmonic as defined in equation 1 can 
be calculated directly using a three-dimensional computer code such as MAFIA. This study uses 
a method where one turn of the helical circuit is simulated. The cold-test characteristics can be 
obtained by applying the quasi-periodic boundary condition of MAFIA to the longitudinal ends. 
This feature of the code permits the user to choose a fixed phase advance per turn in the axial 
direction, allowing the frequency to be obtained at any axial phase shift. The interaction 
impedance is then calculated at this phase shift directly from equation 1 . E^O) is obtained by 

doing a Fourier analysis on the total on-axis axial electric field, £. (0): 
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The time averaged RF power flow P is defined by 


P = wvg , 29 

where v g is the group velocity and w is the time averaged stored electromagnetic energy per unit 
length. 
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with W T the total energy and N the number of helical turns modeled. W x is calculated directly by 
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MAFIA and the group velocity is calculated by taking the slope of the dispersion curve. 


IV. RESULTS 

Using the MMPM helical slow-wave circuit, simulations are compared with experimental 
cold-test measurements below. 

A. Dispersion 

The dispersion was measured experimentally by using the resonant perturbation method. 
An HP 8722D VNA was used to measure the resonant frequencies of the MMPM cold-test 
structure consisting of 159 helical turns [11]. At the helical boundaries are sections of simple 
cylindrical waveguide with diameter equal to that of the slow-wave circuit outer barrel and axial 
length about 1.4 times the length of the longest wavelength of operation. The dispersion was 
also calculated using MAFIA by simulating one helical turn and using quasi-periodic boundary 
conditions. The measured results and those obtained using MAFIA are plotted in Figure 3. The 
agreement is excellent with an absolute average difference across the bandwidth of less than 
0.13%. The CPU time for each frequency point is about 2 hours using a IBM RISC/6000 Model 
595 Workstation. 
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Figure 3 Measured and simulated dispersion for MMPM helical slow-wave circuit 

B. Interaction Impedance 

A cylindrical sapphire rod was placed on the central axis of the MMPM experimental test 
structure and the resonant frequencies measured again. The perturbation rod radius to average 
helix radius ratio is 0.5. The precise positioning of the rod on the center of the structure was 
accomplished by fitting the rod into axial end holes in the test fixture with diameter slightly 
larger than the perturbing rod diameter. Using equations 25 ;and 26, the on-axis interaction 
impedance was calculated and is plotted in Figure 5 along with the results from calculations 
using MAFIA and equation 1 . The results calculated directly from equation 1 using MAFIA are 
consistently lower than measured results with an average dif ference of 26.6%. 

MAFIA was also used to duplicate the perturbation experiment by simulating the helical 
circuit with a cylindrical dielectric rod of size and material p roperties consistent with the 
experimental setup. A cross-section of the MAFIA MMPM helical circuit modeled in the 
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cylindrical coordinate system with the perturbing rod is shown in Figure 4. Computationally, the 
applied frequency necessary to hold the phase shift constant when a perturbing rod is inserted is 
obtained. Using the simulated change in frequency Af, the impedance was calculated in the same 
manner as the experimental values using equations 25 and 26. The MAFIA perturbation 
interaction impedance results are also plotted in Figure 5. These results are coincident with the 
measured values emphasizing the accuracy of the code. Also plotted are the MAFIA 
perturbation results using the first order perturbation expression of equation 27. This data shows 
the large discrepancy between measured results when Lagerstrom’s correction factor is 
incorporated and when it is assumed that the axial electric field is uniform in the transverse 
direction and the transverse fields are neglected. With these assumptions, the interaction 
impedance is an average of 52% lower across the bandwidth compared to measured data using 
Lagerstrom’s correction factor p e . 



Figure 4 MAFIA cross-sectional view of MMPM helical slow-wave circuit with perturbing 

rod 
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Figure 5 On-axis interaction impedance obtained by measurement and using MAFIA with 
the perturbation method (eqn. 25), direct calculation (eqn. 1), and the first order 

perturbation expression (eqn. 27) 


V. INVESTIGATION OF LAGERSTROM’S APPROXIMATIONS 

Figure 5 shows that there is a definite discrepancy between measured impedance and 
results which were obtained using MAFIA and directly calculating the impedance from equation 
1 . This plot, as well as Figure 3, verifies the accuracy of the simulated model by the outstanding 
agreement in the measured and simulated dispersion and imjiedance using the perturbation 
method. The purpose of this section is to investigate the major approximations made in 
Lagerstrom’s analysis and determine their accuracy. To per form a meaningful comparison of 
electric fields for various simulations, it is necessary to compare E 2 /w as each simulation is 
normalized to an arbitrary excitation voltage. For all of the simulations reported on in this 
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report, the energy value stays constant making it valid to compare only the electric fields. It 
should be mentioned that using quasi-periodic boundaries with MAFIA simulates traveling 
waves. Therefore, the fields have amplitude and phase. In the following field comparisons, only 
the amplitudes are shown. 

The major approximations which are made in Lagerstrom’s analysis are: 

A. Only the fundamental axial harmonic of E r is present within the perturbing rod volume (n = 

0 ) 

B. The fields are non-azimuthally varying within the perturbing rod volume (v = 0) 

C. The TM portion of the field dominates the TE portion over the perturbation rod volume (H, 

= E S = 0) 

D. The perturbed and unperturbed electric fields are equal at the perturbing rod surface (E„(b) = 
E,(b)) 

A. Axial space harmonics 

In the analysis of [4] the axial electric fields expressed in equation 7 are assumed to 
contain only the fundamental axial space harmonic, n = 0, over the perturbing rod volume. The 
total longitudinal electric field amplitudes at an azimuthal angle of zero versus axial distance z 
(expressed in number of turns) at radial distances of zero, half the perturbation rod radius (b/2) 
and the perturbation rod radius (b) are plotted in Figure 6 for different values of axial phase 
shift. The wide range of axial phase shifts represents a large frequency range (8.75 - 42.26 GHz) 
so that the reader may get a qualitative idea of the space harmonic content as a function of 
frequency as well as radius. A Fourier analysis was done on the total longitudinal electric field 
along z at the mentioned values of radius and azimuthal angle of zero to give a quantitative 
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analysis of the axial space harmonic content. The amplitudes of the n = 0, ±1 and ±2 space 
harmonics are plotted in Figure 7 (additional space harmonics were not of significant 
amplitude). These plots show that on the axis of the circuit (r = 0), the approximation is 
excellent as the non-zero space harmonic amplitudes are ins ignificant compared to the 
fundamental. However, at finite values of radius the amplitude of the n = -1 space harmonic 
becomes more significant, increasing with radius as one would expect. It is also seen from the 
plots that the amplitude of the n = -1 space harmonic becomes more significant with increasing 
frequency. The inaccuracy of this approximation in the analysis is clearly shown here as the n = 
-1 harmonic is about 50% of the fundamental at the perturbation rod radius at about 42 GHz. 
These results imply that the addition of the n = -1 space harmonic to the measurement analysis 
would be sufficient to improve its accuracy. 



(a) 


20 



BL=90 


25 



z (number of turns) 


(b) 



Figure 6 Total longitudinal electric field amplitude at 0 = 0 as a function of z for (a) (3L=30 
degrees (8.75 GHz) (b) pL=90 degrees (25.72 GHz) and (c) pL=150 degrees 

(42.26 GHz) 
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Figure 7 Magnitude of axial space harmonics of total longitudinal electric field at 0 = 0 for 
(a) PL=30 degrees (8.75 GHz) (b) PL=90 degrees (25.72 GHz) and (c) PL=150 

degrees (42.26 GHz) 


B. Azimuthal space harmonics 

In the analysis of [4] the axial electric fields expressed in equation 7 are assumed to 
contain only the fundamental azimuthal space harmonic, v = 0, over the perturbing rod volume. 
The total axial electric field amplitudes versus azimuthal distance 0 (expressed in degrees) at the 
perturbation rod radius (b) and z = 0.5L are plotted in Figure 8 for several values of axial phase 
shift. From these plots, it is obvious that there is indeed variation in the field with azimuthal 
angle. A Fourier analysis was done on the field along 0 at the perturbation rod radius and z = 
0.5L to give a quantitative analysis of the azimuthal space harmonic content. The amplitudes of 
the v = 0, ±1, ±2 and ±3 space harmonics are plotted in Figure 9 (additional space harmonics 
were not of significant amplitude). These plots show that at the perturbation rod radius, the v = - 
1 space harmonic amplitude is significant equaling about 50% of the fundamental at 42 GHz. 
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The v = 1 space harmonic amplitude begins to emerge at higher frequencies. From inspection of 
Figure 7 and Figure 9, it can be seen that the v' h space harmonic magnitudes are about equal to 
the n ,h space harmonic magnitudes. This verifies the common assumption that v = n because of 
the skew symmetry of the helix [13]. 
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(C) 

Figure 8 Total axial electric field amplitude as a function of 0 at z = 0.5L for (a) pL=30 
degrees (8.75 GHz) (b) |3L=90 degrees (25.72 GHz) and (c) pL=150 degrees 

(42.26 GHz) 
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Figure 9 Magnitude of azimuthal space harmonics of total longitudinal electric field at z = 
0.5L for (a) pL=30 degrees (8.75 GHz) (b) pL=90 degrees (25.72 GHz) and (c) 

pL=150 degrees (42.26 GHz) 

C. TE Fields 

The TE portion of the fields is neglected in the integration of the fields over the 
perturbation rod volume. Thus, from equation 8 the azimuthal component of the electric field is 
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must be zero. The azimuthal electric field amplitudes normalized to the maximum value of 
longitudinal electric field at the perturbation rod radius and 0 = 0 are plotted in Figure 10 versus 
z. From the plot, it can be seen that E e is nontrivial amounting to about 10% of E OTaJ< at PL = 30 
degrees and about 30% at pL = 150 degrees. 
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Figure 10 Azimuthal electric field amplitude normalized to the maximum value of 

longitudinal electric field as a function of z at the perturbation rod radius and 6 

= 0 


D. Perturbed and unperturbed fields 

In the analysis of [4] it is assumed that the perturbed and unperturbed electric fields are 
equal outside the perturbing rod (E„(b) = E,(b)). From this assumption the correction factor p, 
of equation 13 is derived which relates the on-axis longitudinal electric field before and after 

perturbation as E z (0) = p\ E z { 0) . To investigate this approximation, the magnitudes of the 
perturbed fields on the axis obtained directly from MAFIA simulations and as a product of the 
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correction factor p, and the unperturbed field E z (0) are compared. The results are plotted in 
Figure 1 1 showing an average difference of 14% across the bandwidth. 



Figure 11 Perturbed on-axis longitudinal electric field amplitude obtained directly from 
MAFIA simulation and using Lagerstrom’s correction factor p, 

VI. CONCLUSIONS 

The advent of accurate three-dimensional helical circuit models has made it possible, for 
the first time, to investigate standard approximations involved in experimental interaction 
impedance calculations. A summary was given of Lagerstrom’s analysis which relates change in 
axial phase shift or frequency between a perturbed and unperturbed helical circuit to impedance. 
The major approximations were emphasized and followed by an investigation of the accuracy of 
each assumption using MAFIA. The inaccuracy of the approximations along with the large 
discrepancy between measured impedance via perturbation and direct calculations using MAFIA 
prove that a more accurate value of interaction impedance c an be obtained by using three- 
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dimensional computational methods rather than performing costly and time-consuming 
experimental cold-test measurements. This implies a large savings in time and money associated 
with fabricating and testing circuits or scaled circuit models, as well as more accurate results. 

These results are particularly meaningful at higher frequencies where experimental cold- 
testing is more time-consuming, less accurate because of necessarily tight tolerances and more 
difficult to perform. In addition, smaller dimensions associated with higher frequencies imply 
larger b/a values where the approximations in the impedance calculations become least accurate. 
This study also gives an analysis of field behavior in helical circuits which can be used as an aid 
in improving existing helical analysis codes by pointing out which approximations cause 
significant error. 

The accuracy of the MAFIA computational helical model was validated for the MMPM 
helical circuit by comparing the dispersion obtained by measurement and by using MAFIA. The 
agreement was outstanding with less than 0.13% difference. The interaction impedance 
measurements using perturbation and Lagerstrom’s analysis were duplicated using MAFIA and 
compared to measured values with negligible difference. Validation of the computational model 
has also been performed for several other circuit geometries and frequency ranges in [1] 
providing enough confidence in the simulated results to avoid experimental measurements 
altogether. 
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